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ABSTRACT 

We describe the cascade of plasma waves or turbulence injected, presumably by 
reconnection, at scales comparable to the size of a solar flare loop, L ~ 10 9 cm, to scales 
comparable to elementary particle gyro radii, and evaluate their damping at small scales 
by various mechanisms. We show that the classical viscous damping valid on scales 
larger than the collision mean free path (~ 10 8 cm) is unimportant for magnetically 
dominated or low beta plasmas and the primary damping mechanism is the collisionless 
damping by the background particles. We show that the damping rate is proportional 
to the total random momentum density of the particles. For solar flare conditions this 
means that in most flares, except the very large ones, the damping is dominated by 
thermal background electrons. For large flares one requires acceleration of essentially 
all background electrons into a nonthermal distribution so that the accelerated electrons 
can be important in the damping of the waves. In general, damping by thermal protons 
is negligible compared to that of electrons except for quasi-perpendicular propagating 
waves. Damping due to nonthermal protons is also negligible compared to nonthermal 
electrons in most flares which are electron dominated, except for rare proton dominated 
flares with strong nuclear gamma-ray line emission. 

Thus for common flares collisionless damping by background thermal electrons is 
the primary damping mechanism. Using the rate for this process we determine the 
critical scale (or wave vector k c ) below which (i.e. for k > k c ) the damping becomes 
important and the spectrum of the turbulence steepens. This critical scale, however, 
has strong dependence on the angle of propagation of the waves with respect to the 
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background magnetic field direction. The waves can cascade down to very small scales, 
such as the gyro radii of the particles which are well beyond the MHD regime, at small 
angles (quasi-parallel propagation) and possibly near 90 degree (quasi-perpendicular 
propagation). Thus, the spectral distribution would be highly anisotropic at small 
scales. 

Subject headings: plasmas — turbulence — acceleration of particles 



1. INTRODUCTION 

The mechanism of energy release and the process of its transfer to heating and acceleration 
of nonthermal particles in many magnetized astrophysical plasmas in general, and solar flares in 
particular, are still matter of considerable debate. Recent research show that turbulence may 
play an essential role in these processes. In the case of solar flares, it is believed that the energy 
comes from release of stored magnetic energy via reconnection (see Priest & Forbes 2000, Lazarian, 
Vishniac k, Cho 2004 and discussions therein). Both the ordinary and magnetic Reynolds numbers 
Re = l5v/v,Rm = ISv/rj 3> 1, so that the magnetized plasma develops turbulence. Here Sv is the 
velocity change across a turbulent region of scale I, and v and r/ are the viscosity and magnetic 
diffusion coefficients, respectively. More importantly, recent high resolution observations of solar 
flares by Yohkoh and RHESSI satellites have provided ample evidence that, at least from the point of 
view of particle acceleration, plasma turbulence and plasma waves appear to be the most promising 
agent not only for the acceleration mechanism but also the general energizing of flare plasma (see 
e.g. Petrosian & Liu 2004, hereafter PL04, and references cited there). This may also be true in 
other situations (Liu, Petrosian & Melia 2004). These investigations go beyond assuming, e.g. a 
power law electron distribution, as is commonly done, and calculate the expected spectrum based 
on interaction of plasma particles with turbulence. However, in such treatments the spectrum 
W(k) of the turbulence as a function of wavevector k is an input parameter rather than evaluated 
based on first principles. The limitations of such an approach are self-evident. Particle acceleration 
rate depends on the wave spectrum and the wave damping rate are partially determined by the 
particle spectrum. In general, one requires a self-consistent treatment of the coupled wave-particle 
kinetic equations describing the generation of turbulence and its subsequent interactions with the 
background plasma. 

An attempt to solve the problem self-consistently was undertaken by Miller, LaRosa & Moore 
(1996), where coupled equations for energetic particles and turbulence have been studied. For our 
purposes we can rewrite these equations in the following form: 
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Here Dee/E 2 ,A(E)/E and El/E give the diffusion, direct acceleration and energy loss rates of 
the particles, respectively, and Dij(\t)/k 2 and T(k) describe the cascade and damping rates of 
turbulence. The Q's and the terms with the escape times T esc describe the source and leakage of 
particles and waves. Note, that we use a more general form of the equation set compared to that in 
Miller et al. (1996). In our treatment the anisotropy of the turbulent statistics is allowed and the 
leakage of the particles as well as turbulent energy is accounted for. In general, and particularly in 
the case of solar flares that we deal here, it is reasonable to assume that the turbulence is generated 
mainly at a large scale I comparable or somewhat smaller than the spatial extent L of the region 
(with initial velocity SV and magnetic field fluctuation SB). The above equations then determine 
the resultant spectrum and other characteristics of the turbulence, as it cascades to smaller scales 
and is damped by the background thermal plasma, as well as the spectrum of the accelerated 
nonthermal particles. 

In recent years there has been a substantial progress i) in the understanding of cascade of 
incompressible (Goldreich Sz Shridhar 1995) and compressible MHD turbulence (see Cho & Lazarian 
2005a and references therein), and damping of compressible MHD turbulence (Yan & Lazarian 2004, 
henceforth YL04); and ii) in determination of plasma- wave-particle and MHD-turbulence-particle 
interaction rates (see, e.g. Dung & Petrosian 1994; Pryadko & Petrosian 1997, 1998, 1999; PL04; 
Chandran 2000; Yan & Lazarian 2002, 2004, henceforth YL02, YL04, Cho & Lazarian 2005b). 
These advances allow a more thorough description of the coefficients involved in equation (1). In 
this paper we limit our attention to one aspect of this complex problem, namely, to damping of 
turbulence represented by the coefficients T(k). After a brief review of recent progress in the 
understanding of the cascade process in §2 we evaluate the damping rates due to thermal particles 
in the background plasma and a nonthermal population representing the accelerated spectrum 
N(E) (§3). In §4 we summarize our results and discuss their use in our future work on solving the 
coupled wave-particle kinetic equations. 



2. TURBULENCE AND ITS CASCADE 

Plasma turbulence can be decomposed into many wave modes with frequencies lo extending 
essentially from zero frequency to beyond the ion (in our case proton) and electron gyro frequencies 
£l p = eB/rripC and £l e = Q p /S (5 = m e /m p is the electron to proton mass ratio), and wavevectors 
k = l/l spanning the spatial scales from the injection scale to the gyro radius of electrons obeying 
a complex dispersion relation determined by the values of density n, temperature T and magnetic 
field B of the background plasma. 

When we consider turbulence injected at a scale much larger than the proton (or ion) skin 
depth ~ va/^p = 230(10 10 cm~ 3 /n) 1 / 2 cm, where 

(3 A = v A /c = 7 x 10- 3 CB/100G)(10 10 cm- 3 /™) 1/2 (2) 

is the Alfven velocity in unit of speed of light c, initially we deal with modes for which the plasma 
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acts as a single fluid and we are in the MHD regime where the dispersion relation simplifies con- 
siderably. It has been known for decades, that the weak MHD perturbations can be decomposed 
into Alfvenic, slow and fast waves with well-known dispersion relations (see e.g. Sturrock 1994). 
However, it was also believed that such a decomposition does not make much sense for a highly 
non-linear phenomenon of MHD turbulence, where the modes were believed to be strongly coupled 
(see Stone et al. 1999). A study of mode coupling in Cho & Lazarian (2002 and 2003, henceforth 
CL02 and CL03) has shown that the coupling is appreciable only at the injection scale, while along 
the cascade to smaller scales the transfer of energy between the modes is suppressed 3 . This justifies 
the decomposition to different modes even for strong MHD turbulence (see CL02, CL03) and allows 
us to treat their cascade, and interactions with charged particles, of Alfvenic, fast and slow modes 
separately (see YL04, YL03). 

In solar flares, and in many other astrophysical plasmas, one is dealing with a magnetically 
dominated plasma with the plasma beta parameter 

(3 p = MnkT/B 2 = 3.4 x 10- 2 (n/10 10 cm- 3 )(100G/ J B) 2 (T/10 7 K) < 1, (3) 

which also means that the Alfven speed is greater than the sound speed. In this case the slow mode 
can be ignored (see Cho & Lazarian 2005b) and one can use the cold plasma dispersion relation 

lo = VAk cos 9 and u> = VAk (4) 

for the Alfven and fast modes, respectively, where 9 is the angle of propagation of the wave with 
respect to the magnetic field 4 . 

Turbulence generated at large scales can cascade to small scales by nonlinear interactions. 
One important characteristic of turbulence is its self-similarity. Power law spectra were obtained 
numerically for Alfvenic, fast and slow mode turbulence in CL02 and CL03 for the case when 
turbulent energy is injected at large scales. It has also been demonstrated that Alfven (and slow) 
modes exhibit scale-dependent anisotropy similar to that described by Goldreich & Sridhar (1995) 
for incompressible turbulence. 

This can be understood on a qualitative level as follows. For Alfvenic turbulence, the mixing 
motions perpendicular to the magnetic field couple with the wave-like motions parallel to the 
magnetic field providing so-called critical balance condition, k±Vk ~ k^VA, where kit and k± are the 
parallel and perpendicular components of the total wave vector k. This, when combined with the 
Kolmogorov scaling for mixing motions with ~ 5V(k_\_L)~ l l^ (see Lazarian &: Vishniac 1999), 



3 An intuitive insight into this process can be traced to Goldreich- Shridhar (1995) study (see also Lithwick & 
Goldreich 2001, CL02). 

4 In very highly magnetized plasma P\ 3> 1, va in equation (4) should be replaced by va/ yjv\ + c 2 , and the 
actual dispersion relations deviate from these at shorter scales or for k > k crit ~ Q, v / (c0a)(1 + cos 2 6/S) 1 ^ 2 . Fr solar 
flare conditions this scale is larger than the proton gyroradius v t h/Q P ~ 30 cm so that the above dispersion relations 
are good approximations 
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yields a scale dependent anisotropy knL ~ (5V/v A)(k±L) 2 ^ 3 . The mixing motions associated with 
Alfvenic turbulence induce the scale-dependent anisotropy on slow modes, which on their own would 
evolve on substantially longer time scale. The anisotropic spectrum of the Alfvenic turbulence can 
be described as (Cho, Lazarian & Vishniac 2003) 5 

W(k p k ± ) = ^A;- 10/3 exp(-L 1 / 3 |A;|||Af ). (5) 
For Alfven and slow modes the cascade time, same as the hydrodynamic eddy turnover time, 

is 

Teas ^ 1/V k ~ T (k ± L)- 2 / 3 /M A ~ T {k\\L)-\ (6) 

where we have defined a characteristic time and the Alfven Mach number at injection 

t = L/v A = 4.5(L/10 9 cm)(100G/5)(n/10 10 cm- 3 ) 1/2 s and M A = 5V/v A . (7) 

Fast modes in low (3 P plasma, on the other hand, develop on their own, as their phase velocity 
is only marginally affected by the mixing motions induced by Alfven modes. According to CL02, 
fast modes follow an isotropic "acoustic" cascade 6 with W(k) ~ k~ 3 / 2 . 

For such a cascade in each wave- wave collision a small fraction of energy equal to v p h/vk is 
transfered to smaller scales so that the cascade time scale is characterized by (CL02): 

Teas = (V ph /V k )(l/V k ) = T^kLy 1 ! 2 / M\. (8) 

Here v p h = oj/k = va is the phase velocity of the fast mode, and we have used the scaling relation 
Vk = 5V(kL)~ l l A appropriate for fast modes. 

For solar flare conditions f3 p < 0.1, and /3a = va/c ~ 1CT 2 . Assuming Ma < 1 the above 
cascade times are about few seconds at the injection scale but much shorter at shorter scales. 

Alfven and slow modes are inefficient in scattering energetic particles (Chandran 2000, YL02). 
YL02 identified isotropic fast modes as the dominant scattering agent. It is also possible to show 
(see YL03) that fast modes are the dominant mechanism for acceleration of particles via resonant 
interaction. Fast modes also dominate the acceleration of particles through non-resonant interaction 
(Cho & Lazarian 2005b). Consequently in what follows we will concentrate on the damping of the 
fast modes. 



5 Note that integrating over the parallel and perpendicular components one gets W(k±) oc fc ± J ^ 3 and W(k^ ) oc fey 2 , 
respectively. 

6 This scaling, plausible on theoretical grounds (CL03, Cho & Lazarian 2005a), is supported by numerical simula- 
tions of CL03, but may still require more testing. On the other hand, in the high (3 P case the scaling of the fast modes 
trivially reduces to the scaling of acoustic turbulence because the fast modes are essentially acoustic fluctuations. 
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3. DAMPING RATE OF TURBULENCE 

The second important process determining the spectrum of turbulence is its damping rate. 
Damping becomes important whenever the damping time T~ 1 (k) becomes comparable to or shorter 
than the cascading time T cas (k). As we shall show below, for solar flare conditions, the damping 
time is longer than the cascade time at large scales but decreases faster with decreasing scale and 
becomes dominant above the critical wave vector where T(k c )T ca s(kc) = 1- In this section we derive 
the damping rate and the critical wavenumber k c . We first describe the damping by the background 
thermal plasma. 



3.1. Thermal Damping 

In fully ionized plasma, the damping can be divided into two parts: collisional and collisionless 
with their regimes of relevance determined by the ratio of the turbulence scale and the Coulomb 
collision mean free path of the background plasma (Braginskii 1965), 

7 / T \ 2 /10 10 cm- 3 \ 

Viscous damping is important for scales I > Acoub so that it can play a role between the injection 
scale < L ~ 10 9 cm and Acoul ~ 10 8 cm. For smaller scales, fcAcoul > 1, the damping rate is 
determined by less efficient collisionless processes. 



3.1.1. Viscous Damping 

The viscous damping rate is derived in Appendix A where we show that for low beta plasma 
of interest here we have 

T ms (k,e) = 0.13T -Vp /2 (^)(kA C oui)sin 2 fl, for k\ Co ui < 1. (10) 

By equating r~ 1 (/c) from above equation to T cas in equation (8) we obtain the critical scale or 
wavevector 

Mcou! = 3.9(M A /sin#) 4 / 3 (A Co ui/A3 P ) 1/3 - (11) 
For Ma ~ 1 and (3 P < 0.1 the last two terms are greater than one indicating that the critical 
scale is less than the Coulomb mean free path where this damping rate is not valid; k c \cou\ < 1 if 
sin# > 2.8Ma, so that viscous damping could be marginally important for Ma < 0.3. 



3.1.2. Collisionless Thermal Damping 



The nature of collisionless damping is closely related to the radiation of charged particles in 
magnetic field. Charged particles can emit plasma waves through acceleration (cyclotron radiation) 
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and Cerenkov effect, and can also absorb the radiation under the same condition and cause damping 
of the waves (Ginsburg 1961). For example the gyroresonance with thermal ions causes the damping 
of the modes with frequencies close to the ion-cyclotron frequency (Leamon, Mathaeus & Smith 
1998). The particles can also be accelerated either by the parallel electric field (Landau damping) 
or the magnetic mirror associated with the comoving compressible modes under the Cerenkov 
condition k\\v\\ ~ u, known as transit time damping, or TTD for short. Because head-on collisions 
are more frequent than trailing ones, the energy is transfered from waves to particles. 

For small amplitude waves, particles should have parallel speed comparable to wave phase 
velocity to be trapped in the moving mirrors. This gives rise to the above Chernenko condition. 
For a thermal plasma, this requires thermal speed v t h ~ v P h or a plasma beta of oder unity. At 
lower values of (3 P the fraction of particles satisfying this condition is smaller which means that 
the damping rate decreases with decreasing (5 p . The damping rate of fast mode with frequency 
ijj = kvA and f3 p < 1 can be written as 



T th (k,9) = r x 



exp I — — — —5— I + -^= exp 



9(0), (12) 



/3 p cos 2 9j y/s V /3 P cos 2 9 / 

where, as stated above 5 = m e /m p is the electron to proton ratio and we have defined a characteristic 
damping rate 

_ y^pp~5 sin 2 9 

T = (kL) -. (13) 

2to cos 9 



This damping rate, without the last term g{9) and valid for 9 > \Ju/VL p , coincides with the one in 
Ginsburg (1961). 7 In the square bracket the first term represents the contribution from electrons, 
and the second term is due to protons. The function g{9) for 9 <C 1 is 



m 1 i + - t ,('■ «i»">^A¥; (14) 



2 



derived by Stepanov (1958) extends the relation to small angles where the damping rate decreases 
by a factor of two. (In the limit of 9 = 0°, there are no compressions, fast modes are degenerate 
with Alfven modes.) For (3 P < 0.1 and sufficiently large 9, the damping due to electrons dominates 
and the damping rate can be written in a simple form 

T th (k,9) =r exp[-<V(/3 P cos 2 #)], for fcA Co ui > 1, (15) 

where we have ignored the correction g{9) at small angles. Note the similarity of this relation with 
that for viscous damping; the main difference is the absence of the extra term (fcAcoul) in equation 



Note that this expression differs from the one in Ginsburg by a factor of 2. This is because we are concerned 
with the energy instead of wave amplitude damping. 
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(10). One can then combine the two expressions to obtain an approximate damping rate valid at 
all scales 



where we have deleted the exponential part in equation (15) which is equal to one except at small 
range of angles near 7r/2 (i.e. cos 9 < .023/\//L). A more accurate expression is obtained if one 
divides Q by the square bracketed term and g{9) in equation (12). 

In Figure (1) we compare the cascading time with the damping time (r^ = 1/Tth) at dif- 
ferent scales for (3 P = 0.01 and different angles (right panel), and for 6 = 45° at three values of 
(3 P = 0.001,0.01 and 0.1 (left panel), corresponding to magnetic fields of B ~ 600,180 and 60 
G, respectively (see eq[3]). We use typical solar flare values; temperature T = 10 7 K, density 
n = 10 10 cm~ 3 , and we set Ma = (SV/va) = 0.3. The angular dependence enters this damping 
by two competing factors. In general the damping increase with because magnetic compression 
increases so that more particles can be trapped and interact with the waves. However, when 8 
approaches to 90°, i.e. for quasi-perpendicular propagation, most thermal particles will not be in 
resonance with the fast mode waves in a low [3 P medium, which explains the decrease of damping 
in this regime. 

By equating the collisionless damping time from equation (15) with the cascade time in equa- 
tion (8), we attain the critical wave vector 



The variation of k c with angle for the thermal collisionless damping using the exact expressions 
is shown in Figure (2). As evident the damping scale given by equation (17) varies considerably 
especially when 9 — ► and 9 — ► 90°, where it becomes smaller than the proton gyro radius (shown 
by the dashed line). Note also that for (5 P < 0.1 the damping scale is larger than the collision mean 
free path (or fc c Acoui > 1, shown by the dashed-dot horizontal line), except for few degree around 
9 = 85°, which is within the range of its validity. The specific range of the 9 where the relations 
breakdown depends on the plasma (3 P . 

This describes the well known fact that at large k vectors (i.e. small scales) the turbulence 
will be very anisotropic. Only quasi-parallel and quasi-perpendicular modes survive at such large 
wavevectors (and corresponding frequencies) which are needed in the acceleration of low energy 
particles. However, as turbulence undergoes cascade and/or waves propagate in a turbulent medium 
the character of this anisotropy changes because the angle 9 is changing due to the randomization 
of wave vector k and the wandering of the magnetic field lines discussed next. 



r tot (k,9)~r., 



WMVCl+C 1 ), ^ C = 6.3cosfl(fcA C oui), 



(16) 




(17) 
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Fig. 1. — Variation with wavevector k (in units of injection scale L" 1 ) of the cascade and damping 
timescales (in units of tq = L/va) of fast modes at 6 = 45° for different values of f3 p (left panel) 
and for various 6 for p = 0.01 (right panel). On the left panel the short dashed, long dashed and 
solid lines are for viscous damping, collisionless damping and the total damping, respectively. Note 
the steepening of the total damping at low values kL < L/Acoul (represented by the left vertical 
dotted line) due to emergence of viscous damping. The critical damping wavevector k c is given by 
the intersection of the two time scales, which as shown in Figure (2) depends on the angle 9. The 
vertical dotted lines at the right side show the scale of the proton thermal gyro radius L/r g ^ p 



-10- 



10 



15 



Damping lengthscale vs. 



10 L 



(3=0.01 , collisionless 
p=0.1, collisionless 
mean free path 
proton gyro scale 




10 20 30 40 50 60 70 80 90 

9 



Fig. 2. — Normalized cutoff scales of turbulence due to thermal damping vs. the angle 9 between 
k and B obtained by equating turbulence cascading and damping rates from equations (8) and 
(12). The horizontal dash-dotted and dashed lines represents the scales of the mean free path 
(L/Acoul)) the dividing line between collisional and collision less damping, and the thermal proton 
gyro radius (L/r gtP ), the limit of applicability of MHD, respectively. The dotted lines show the effect 
of field wandering and k vector randomization at angles near 90° for (3 P = 0.1. The wanderings are 
ineffective at (3 P < 0.05. Note that the critical scale is almost always smaller than the mean free 
path so that the viscous damping can play only a marginal role (see text). 
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3.1.3. Damping Anisotropy 

The damping of fast modes described above is valid for small perturbations and for a uniform 
background magnetic field. A more realistic setting for damping in turbulent media, which is based 
on better understanding of turbulent cascade (Cho & Lazarian 2002, 2003) and magnetic field 
wandering (Lazarian & Vishniac 1999, Lazarian, Vishniac & Cho 2004), is discussed in YL04 for 
the case of interstellar medium. Here we apply this approach for turbulence in solar flares. 

For fast mode cascade, the non-linear cascading occurs by interaction of wave packets that 
are collinear (see review by Cho &; Lazarian 2004), to = lo\ + 002 and k = ki + k2, and the final 
wave is parallel to the initial ones. However, the energy conservation is true within an uncertainty 
condition 5lo ~ t^J, <C uj, which will give rise to an orthogonal wavevector component. It is easy 
to show then that this uncertainty in frequency will yield a small transverse component 5k related 
to 8u> as 800 ~ V p h8k(Sk/k). For Slo/oj ~ \/(uoT cas ) <C 1 the angle 50 between the final and initial 
wavevectors will be small. Combining this with equation (3) we get for a given Alfvenic Mach 
number Ma the following expression (YL04) 

tan 56 ~ 5k /k ~ {5uj/uj) 1 / 2 ~ M A (kL)~ l / A , (18) 

It is easy to see that, as in interestellar medium case, the field line wandering in the case 
of solar flares is mainly caused by the shearing due to Alfven modes. According to Lazarian & 
Vishniac (1999) 8 the field line diffusion along and perpendicular the mean field provides 

tan«50|| ~Mi(— J , tan<50 ± ~ M A (-J (19) 

where z and r are the distances along and perpendicular to the mean field direction. During one 
cascading time (equation 8), the fast modes propagate a distance r cas = T^j [M\\fkV) and see an 
angular deviation 

-«*>/-°».+-'«W#iW^f <20> 

Note that at the largest scale kL = 1 the randomization is of the order of Ma < 1 9 and 
decreases slowly with k; 59 scales as (fcL) -1 / 4 at small angles or quasi-parallel modes and as 
~ (/cL) -1 / 6 at other angles. This shows that randomization decreases with the decrease of the 
scale. 



8 A typo of numerical coefficient has been corrected. 

9 It is useful, in general, to define the domain of MHD turbulence as starting from Alfven Mach number of unity 
(see discussion in Cho & Lazarian 2005a) . This condition is automatically satisfied for strongly magnetized turbulence 
in solar flares. 
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While for the processes of scattering and acceleration of particles by fast modes randomization 
is important (Yan & Lazarian 2004), it does not always result in tangible changes of the overall 
picture of damping. For instance, combing equations (17) and (20) we can show that for quasi- 
parallel modes (i.e. 9 ~ 0) at the critical wavevectors k c the randomization angle is rather small, 
i.e. 59 ~ 10 -3 Ma$J 2 . For other angles 59 < (/? p /30) 2 / 3 and is still small but may not be negligible. 
In particular, for the quasi-perpendicular modes, we can make an estimate of the presence of field 
line wandering on the damping truncation scale by evaluating the average of the damping rate in 
equation (12) over the small range tt/2 — 59 to vr/2, where from equation (20), 59 ~ M^ 3 (feL) -1 / 6 . 
For 59 <C 1 we can define a = ir/2 — 9 and use the approximations sin 9 = 1 and cos 9 = a so that 
the average value of the damping rate given in equation (12) is roughly given by 

V^p(fcL) [*> y^ggfcL) / 5 \ 

{Tth) = 2r 59 J, da 6XP - 5/iPpa ),a = 4r 59 El ) ' (21) 

Equating this with r cas from equation (8) we get x 2 Ei(x) = A with x = (5 / f3 p )(kL / M\) l l 3 and 
A = (4/^)M A l 5 3 / 2 Pp 5/2 = (.025//3p)~ 5 / 2 , where we set M A = 0.3. For (3 P = 0.1 this relation is 
satisfied for x ~ 0.15 or k c L = M^(x/3 p /5) 3 ~ 20. This scale is shown by the horizontal dotted 
line in Figure (2). This means that modes in the cone near 90° get damped above this scale due to 
randomization of 59. However, because x 2 Ei(x) < 0.22 there is no solution for f3 p < 0.05 indicating 
that the field wanderings do not reduce the damping scale at 90° for highly magnetized plasmas. 

These results are only rough estimates. A more detailed study of the issue is more appropriate 
in the context of particle acceleration in solar flares and will be done elsewhere. 



3.2. Nonthermal damping 

The wave damping rate calculated above assumes that the energy lost by waves with spectrum 
W(k) goes into heating the plasma and that the plasma maintains its Maxwellian distribution, 
presumably via Coulomb collisions in a timescale of tc ou i ~ ^Coui/cs, where cs is the sound 
speed. This requires a longer damping time; T^ 1 > rc ui- As stated above we are interested in a 
collisionless plasma with r cas <C 7c ui- This combined with the fact that damping is important when 
T^ 1 < T cas implies that the above condition is not satisfied and some particles get accelerated to 
energies much higher than ksT. Because Tcoui decreases with energy fairly rapidly it can be shorter 
that the other times at low energies, the particle spectrum there will be approximately Maxwellian. 
As shown in PL04 solution of the particle kinetic equation (1) with a given background thermal 
plasma and an assumed spectrum of turbulence does lead to a particle spectrum consisting of a 
quasi-thermal part and a nonthermal tail with dividing energy roughly where Tcoui is equal to the 
acceleration time r ac ~ E/A(E). Thus, we need to also consider damping of the waves by the 
nonthermal tail. As mentioned at the outset one must carry out this self-consistently by solving 
the coupled wave-particle kinetic equations which is beyond the scope of this paper. Here we derive 
the damping rate due to electrons and protons with a total density of iVo, a power law energy 
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spectrum, N(E) = N (a — l)(E/E ) a /E for E > E , and isotropic pitch angle distribution. 

As can be surmised the calculations of the particle diffusion coefficient D pp = ( A ^ p ), or 
the acceleration rate A(E)/E and the damping rate r nont / l (k) are intimately connected. Let us 
represent the transition rate (integrated over the particle pitch angle cosine fi) of the interaction 
between a wave k and a particle with energy E by <r(k, E). The Fokker-Planck diffusion coefficients 
D(p), which is D pp integrated over fi, is 

dfiD pp = / d 3 kW{k)a(k,E). (22) 

-1 J k m i n 

From this we can get the rate of systematic energy gain by particles (PL04): 

M^-^^-T (23) 

which also means that 

Because the energy lost by turbulence \V nc ,nth = / InonthityW (k)d 3 k to nonthermal particles must 
be equal to the energy gain by these particles, £ = J A(E)N(E)dE, the damping rate is given by 

T nonth (k)= dEN(E)X(k,E). (25) 
J Eo 

The transition rate a for a particle with gyrofrquency Q and pitch angle cosine \x interacting 
with a wave of frequency uj and wavevector k is determined by the resonant condition 

u-kcos6vn = nn/-/, n = (0, ±1, ±2, ...). (26) 

This rate can be expressed as sum over n of squares of Bessel functions J?(kj_vj_ , y/£l) with i = 
n, n — lorn + 1 and v± = v(l — fi 2 ) 1 ^ 2 . (For details see Pryadko & Petrosian 1999.) For waves 
propagating parallel to the B field k± = and only Jo(0) / and n = ±1 (n = term also 
vanishes) . The resonance condition then requires k res ~ Qjvn ~ t„ . Given the parameters we 
adopt here r g ~ 1 and 50 cm for thermal electrons and protons, respectively, which is certainly 
beyond the MHD regime. Only quasi-parallel propagating waves cascade to such small scales 
without undergoing thermal damping and can contribute to acceleration of low energy particles 
(see Fig. 2) 

For obliquely propagating waves (fast or Alfven) things are more complicated and all n's could 
contribute. However, these waves are damped at scales much larger than that required for the above 
resonant condition. Except for n = or the transit time damping (or TTD for short) mode which 
happens at all scales with the resonant condition v/u, = va/ cos 8. TTD is resonant interaction with 
parallel magnetic mirror force. Thus, in what follows we consider this process with the transition 
rate 
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a(KE) = ^- f ^(1-M 2 )f -, ( T r l (27) 

2nm p J_ 1 y Tca i + {k\\vn-oj) 2 J ft 

where n is the density of protons (which we take to be equal to that of electrons) , and The resonance 
function in the large parenthesis is produced by integration over time. In general, the width of this 
function Afi = T~ a l/k^v = M\(kL)' 1 / 2 (v a/ v cos 9) <C 1 because usually Alfven Mach number 
Ma < 1 and because damping normally becomes important at scale k^ 1 <C L. (This is not true 
only for nearly perpendicular propagation.) Then the resonance function can be approximated 
by a delta function, i.e., t~^ s I{t~ 2 s + (k\\v/j, — uj) 2 ) — > nS^kuVfj, — u) = tt5((i — (i res )/k\\v, where 
[i res = va/v cos9 . In Figure (3), the integrand of cr(k, E) is plotted versus jjl for some interesting 
cases. As evident, when the resonance condition is satisfied, i.e. v\\ < va cos 9, the transition rate a 
peaks sharply by many orders of magnitude so that it can be well represented by a delta function. 

With this simplification the integration over [i can then be carried out easily. Then from the 
relation J[ = (Jq — J2V2 and equations (27), (25) and (24) we can obtain the damping rate 

7r Q 2 m f R 2 \ f°° 

r nontfc (k) = 1 - % / dEN(E)Q(E - E c )/((3 7 ) 

8 nnipckr] \ V J Je 

[2jf(x)+xJ 1 (x)( j (x)-J 2 (x) )}, x = p-yckSl-Wl ~V 2 ^~ Pa/(Pv) 2 - (28) 

Here we have defined rj = cos 9 and E c /mc 2 = (l/yjl — P A /rj 2 ) — 1, and @(x) = 1 for x > 1, and 
zero otherwise, is the Heaviside step function. 

Damping by Electrons: For electrons kc/Q e ~ 10 _8 (/cL)(100G/i?) <C 1 for the relevant 
scales k < k c . Therefore, except for extreme relativistic electrons the variable x <C 1 and we can 
use the first order approximation J n (x) ~ (x/2) n /nl for the Bessel functions. Then the terms in 
the square brackets in the above equation can be approximated as x 2 /2 to give 

where E m = m&x(Eo, E c ). 

We first note that if we use a nonrelativisitic Maxwellian distribution, N(E) = n(2/ y/Tt)(kBT)~^/ 2 E 1 / 2 exp (— E 
this gives the damping rate of 

r.,« = r„exp(-^) (30) 

which is identical to the electron part of the collisionless thermal damping given in equation (12). 

For a nonthermal distribution we can carryout the integration which leads to a complicated 
expression shown in Appendix B. In the nonrelativistic and extreme relativistic cases the result 
simplifies considerably. For the spectral index a > 2 most of the contribution to the integral comes 
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Fig. 3. — The dependence of the diffusion coefficient on \x from equation 27 for 6 = 80° (left 
panels) and 6 = 88° (right panels), and for particle kinetic energies of lOkeV (top panels) and 
5MeV (bottom panels). Most resonance can be approximated by 5 function (see text). In those 
cases -ucos^ < va, the resonance condition is not satisfied and therefore there is relatively very 
weak interaction rate. 



from the low energies so that when E m <C m e c 2 (i.e. Eq <C m e c 2 and E c /m e c 2 ~ P\/(2r] 2 ) <C 1) 
we can use the nonrelativistic approximation to obtain 



(31) 

The second expression is valid near 9 = n/2 where r/ <C Pa/s/Eo- For the solar flare conditions 
we are considering here, Pa ~ 0.007 and Eq ~ 0.02m e c 2 ~ 10 keV, this expression is applicable 
only for tt/2 — 9 < 0.05 or within 0.3° of perpendicular direction. Thus, for all angles outside this 
range we have the first expression which is different from the thermal damping rate of equation 
(30) by the presence of the terms in the square brackets. The main part here is the ratio Nq/u 
of the nonthermal to thermal particle densities and by (-Eo/^bT) 1 / 2 , i.e. the ratio of the mean 
momenta (or velocities) of nonthermal to thermal electrons. In general, this term is less than one 
and damping by thermal electrons is dominant. However, in some large solar flares one requires 
acceleration of a large fraction of the background thermal particles to energies 3> ksT so that the 
damping by nonthermal particles could be significant. 

If Eq 3> m e c 2 , we can use the extreme relativistic approximation to obtain 

^-.( 1 -f) ! ©" 2 (?)^)^)" , (^)"' (32 » 

which is valid for all angles, where a E c /m e c 2 ~ Pa/v 2 ^ EQ/m e c 2 , so that E m = Eq » m e c 2 , 
except for a extremely narrow range of angles give by (1 — (m e c 2 / Eq) 2 P A Ci| 2 < P\. One can 
combine the above two expressions to obtain an approximate relation valid at all energies: 

r„ mi ak) = r = '(i + ^|(^) V y (33) 

However, it should be noted that the extreme relativistic equation is valid only for scales k < 
/ [C70 sin 9^1 — P A /r] 2 ], otherwise the approximation used for the Bessel functions breaks down. 
For most angles, and for parameter values adapted here, this means kl < 5 x 10 7 /7o- But for 9^0 
or cos 9 — ► Pa the above expression would be valid at much smaller scales or larger values of kL. 
These limitations are also true for the more general equation (Bl) in the appendix. 

Damping by ions: For protons (and other heavier ions) the condition k±v±/Qi <C 1 is 
not always satisfied. Nevertheless, if we use the small argument asymptotic expression J n (x) ~ 
(x/2) n /n\ for Bessel functions (as done above for electrons), we can get similar estimate for the 
damping rate due to interaction with protons. 

For example, for a Maxwellian proton distribution one can show that the resultant damping 
rate will be same as that for electrons with m e — > m p , which means setting 5 = 1 in the equation 



jvnonrel 
1 nonth 



(k) 



1 - 



Pa 



( ^(a-\)No \ f Eq \ 
V 4(a-3/2)n J \k B T ) 



1/2 



-17- 



(30). Aside from the factor of 5 this is identical to the contribution by protons to the collisionless 
thermal damping in equation (12). 

Within a similar accuracy we can also estimate the damping rate due to nonthermal protons. 
Ignoring angles near ir/2 for the moment, from equations (28) or (29) we note that T non th oc (p)Nq 
so the relative importance of protons and electrons depend on their total number ratios and their 
mean momenta which will be same as the ratio of the momenta po at the low end of the spectrum. 
In solar flares we deal with nonrelativistic values of Eq = y/2mpo and much fewer number of 
accelerated protons compared to electrons. In majority of flares the ratio of the total energies 
R = (N E ) p /(N E ) e (in the observable range (E 0t€ > 10 keV for electrons and E QtP > 10 MeV 
for protons) is much less than one and varies from 0.01 to 10 in flares with detectable gamma- 
ray line emission produced by the accelerated protons and ions (Miller et al. 1997, attributed 
to R. Ramaty & N. Mandzavidze) . Thus T nont h,p/T nont h,e = R\/ -^o.e/ (Eo tP 5) ~ R for the above 
mentioned energies. In summary, usually one can ignore the nonthermal damping due to protons 
relative to electrons (as was the case for the collisionless thermal damping) in most solar flares. 
But nonthermal damping by protons could be more important than that of nonthermal electrons 
in flares with strong gamma-ray line emission. 

4. SUMMARY AND DISCUSSION 

The study of interactions of plasma waves and turbulence with particles of magnetized plasma 
is a complex process and requires an accurate formulation of the cascade of the turbulence from 
a large injection scale to smaller scales and the damping of the waves by the background thermal 
particles and those accelerated into a super-thermal power-law tail arising from these interactions. 
Equipped with this knowledge one can then determine the evolution of the spectrum and angular 
characteristics of the turbulence and particles by solving the coupled kinetic equations (see eq.[l]). 
There has been considerable progress in the understanding of the cascade of turbulence from an 
injected scale L to lower scales or higher wavevector k and its expected spectral and angular 
characteristics. We briefly review these for magnetically dominated or low beta plasmas ({3 P < 1), 
such as those envisioned for solar flares, and indicate that of the three MHD modes Alfven, slow 
and fast, the latter can play a dominant role in heating and acceleration of plasma particles. The 
aspect of this process most relevant to our goals in this paper is the rate of cascade or cascade time 
as a function of scale of the turbulence. In general, the cascade time T cas for all these modes is of 
order of Alfven injection scale time To = L/va and decreases as (kL)^ 1 ^ 2 . 

The spectral and angular distribution of the turbulence is further modified at smaller scales 
when the damping rate becomes comparable and larger than the cascade rate. The main goal of 
this paper is to give a complete description of the damping process. We review the basic processes 
involved here and present equations describing the damping rate of the turbulence due to different 
mechanisms. We first consider viscous damping valid on scales larger than the collision mean 
free path, or for /cAcoul < 1, with the damping time scale ~ To[\f]3 p (kL)(k\Q OVL \)\~ l oc k~~ 2 . 
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Because of this rapid decline this damping can become quickly important and stop the cascade 
process. This would be the case for high beta plasmas but for solar flare conditions, this damping 
mechanism is applicable in the narrow range of scales between the injection scale L ~ 10 9 cm and 
Acoul ~ 10 8 cm, where because of small value of f3 p , > r C(ls and can be neglected. For smaller 
scales the damping is produced by collisionless processes. Here we have described the damping due 
to thermal (Maxwellian) and a nonthermal (power-law) distributions, separately. 

The thermal damping is dominated by electrons and for most practical purposes can be ap- 
proximated as Td ~ To[y/]3 p {kL)] _1 oc k^ 1 . Proton contribution to this process can be important for 
PpCos 2 6 > 0.18 which will not be the case for low beta plasmas (3 P < 0.1 under consideration here. 
Here 9 is the angle between the magnetic field and the propagation direction. The same is true 
for other ions. We combine the collisionless damping valid for fcAcoul > 1 with that for the viscous 
damping and give a simple expression valid approximately at all scales. Equating the damping and 
cascade times we determine the critical wavevector k c above which the damping becomes dominant 
and would cause the spectrum of the turbulence to steepen. The damping is highly anisotropic and 
the critical wavevector varies considerably with angle 6, being much larger for quasi-parallel and 
quasi-perpendicular propagations. However, we show that this anisotropy around the perpendicular 
direction is smoothed out by magnetic field wanderings caused by the shearing due to Alfven modes 
for f3 p < 0.05. The quasi-parallel waves, on the other hand, are not affected by this process and 
can survive without damping to scales as small as the particle gyro radii where the MHD regime 
breaks down and other plasma and kinetic effects become important. These parallel propagating 
waves may be then the most important agents for acceleration of low energy electrons, protons and 
other ions (see PL04, and Liu, Petrosian & Mason 2005a, 2005b). 

We have also evaluated the collisionless damping rate due to a population of nonthermal elec- 
trons and protons. We have argued that the most important process here is the transit time damping 
mechanism, and show that this process gives a damping rate very similar to that obtained for a 
thermal distribution. In general, the damping rate is essentially proportional to the mean momen- 
tum times the number of the particles. Thus the relative importance of thermal and nonthermal 
populations depends on the product of the ratios of their densities and average momenta. In most 
cases except for extremely hard nonthermal tails (electron index a > —1.5 or —2, for nonrelativis- 
tic and extreme relativistic cases, respectively), this ratio will be less than the energy content of 
the two population. In particular, this will be true for most solar flares except for the strongest 
bursts where one requires acceleration of all the available background electrons. This behavior also 
indicates, that as is the case of thermal damping, here also the contribution of protons relative to 
electrons can be neglected except for very rare flares with strong nuclear gamma-ray line emission 
which require more energy for accelerated protons than electrons. Most flares, however, are elec- 
tron dominated and the contribution of nonthermal electron will increase the damping rate by the 
above basic ratio at all k and 9 and decrease the value of the critical wavevector but not affect its 
anisotropy. 

We stated that fast modes dominate slow and Alfvenic modes in terms of acceleration and 
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therefore concentrated on the fast modes. This is only true if turbulence is strong, i.e. when the 
critical balanced condition is satisfied for Alfvenic modes. At the injection scale the turbulence may 
be weak (see Galtier et al. 2000) and develop a cascade with ku = const. However, such a cascade 
has a limited inertial range (see discussion in Lazarian & Vishniac 1999, Cho, Lazarian & Vishniac 
2003) and beyond this range it transfers to the strong Alfvenic turbulence. Therefore disregarding 
weak turbulence seems justified 10 . Moreover, reconnection events should produce perturbations 
with velocity of the order of Alfven velocity, which should produce strong turbulence from the very 
beginning. 

MHD turbulence that we considered was balanced in the sense that the equal flux of energy was 
assumed in every possible direction. In solar corona we expect the energy injection to be localized 
both in space and in time. As the result turbulent energy propagates from such sources, e.g. 
reconnection region creating an imbalanced cascade. The properties of imbalanced turbulence (see 
Maron & Goldreich 2001, Cho, Lazarian & Vishniac 2002, Lithwick & Goldreich 2003), in particular 
its damping time and scaling, can be very different from the balanced one. The Alfvenic cascade is 
being strongly modified by imbalance and this may result in much lower rates of cascading, if the 
imbalance is strong. However, variations in Alfven speed that are present in solar corona are likely 
to result in reflecting Alfven perturbations. These reflections mitigate the imbalance and therefore 
we believe that the effects of imbalance will not be substantial. A more detailed study of the issue 
will be presented elsewhere. 

We have limited our considerations to scales above the particle gyro radii, where the MHD 
approximation is valid. For shorter scales one must consider mechanisms of damping or acceleration 
other than the TTD. In particular gyro-resonance scattering must be included with more realistic 
dispersion relations than those given in equation (4). We intend to address these extensions of the 
current results in our future publications. 

Transfer of turbulent energy from large to small scales and its damping is a general process 
that can be important for heating and particle acceleration in various environments other than 
solar flares, such as in gamma-ray bursts (see Lazarian et al. 2003) and accretion around black 
holes (Liu, Petrosian & Melia 2005). 

We acknowledge support from the NSF grants ATM-0312282 at Wisconsin and ATM-0312344 
at Stanford. AL also acknowledge the NSF Center for Studies of Magnetic Self-Organization in 
Laboratory and Astrophysical Plasmas, AL and HY acknowledge partial support from NSF AST- 
0307869, and VP would like to acknowledge support from NASA grants NAG5-12111 and NAG5- 
11918-1. 



1 More definite statements can be obtained if velocity fluctuations associated solar flares are analysed. The cor- 
responding techniques developed for the interstellar medium, e.g. Velocity Chanel Analysis, Velocity Correlation 
Spectrum (Lazarian & Pogosyan 2000, 2004), Modified Velocity Centroids (Lazarian & Esquivel 2003, Esquivel & 
Lazarian 2005) can be applied for the purpose. 
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A. Damping Due To Ion Viscosity 

Viscous damping is not important unless there is compression. Therefore it only influences 
compressible modes and has marginal effects on Alfven modes. In a strong magnetic field, the proton 
gyrofrequency is much larger than proton collisional frequency j ~ c^/Acoub (^p r Coui = 5 x 
10 3 (i3/100G)(10 10 cm~ 3 /n)(T/10 7 K) 3//2 ), and the transport of transverse momentum is prohibited 
by the magnetic field. Thus, transverse viscosity coefficient rj± ~ ?7o/(^pi"Coui) 2 is much smaller 
than longitudinal viscosity coefficient t]q = 0.96n/csTTc O ui (see, e.g. Braginskii 1965). 

Considering only the zeroth order terms due to longitudinal viscosity, the viscosity tensors are 
k xx = nyy = -T] (W XX + W yy )/2 and tt zz = -r] W zz , where W jk = ^ + |^ - §5,- fc V . v is the 
rate-of-strain tensor, v is the fluid velocity (Sigmar 2002). Here, z axis is defined by the magnetic 
field. Heat generated by the viscosity is 



du dv dv 

Q vis = tt : Vv = -tt xx -^ - Kxx-q 1 ~ ^zz-q^ = Vo(dv x / dx + dvy/dy - 2dv z /dz) 2 /3. (Al) 

Dividing this by the total energy associated with the fast modes, we obtain the damping rate T V i S . 
While the damping due to compression along the magnetic fields (the 3rd term) can be easily 
understood, it is somewhat counterintuitive that the compression perpendicular to magnetic field 
also results in damping through longitudinal viscosity. However, the origin of this viscosity can 
be easily traced (see Braginskii 1965). Indeed, for motions perpendicular to the magnetic field B, 
• v = hjn ~ B/B, implies the transverse energy of the ions increases due to the adiabatic 
invariant v\jB. If the rate of compression is faster than that of collisions, the ion distribution 
in the momentum space will become distorted away from the isotropic Maxwellian sphere to an 
oblate spheroid with the long axis perpendicular to the magnetic field. As a result, the trans- 
verse pressure becomes greater than the longitudinal pressure by a factor Tcouin/n, resulting in 
a stress ~ PTQ ou \h/n ~ r/oVj_ • v, where P = nksT is the longitudinal pressure. The restora- 
tion of the equilibrium increases the entropy and causes the dissipation of energy. In a low (5 P 
medium, compressions are perpendicular to magnetic field, thus T v i s = k^rjo/Snrrii. In high (3 p 
medium, as pointed out in §2, the velocity perturbations are radial. Thus according to Eq.(Al), 
the corresponding damping rate T v i s = /c 2 ??o(l — 3 cos 2 0) 2 /(Snrrii). 

Putting all these together for rrii = m p we get 

or in terms of physical parameters appropriate for solar flares 

r OOSc-V/^P^V T V7 l0 10 cm- 3 ysin 2 g for (3 p « 1 

V ms - 0.05s (kL) ^ L j \ v io 7 K / V n Jl(l-3cos 2 e) 2 for p » 1 (A3) 
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This damping becomes important at a scale where T v i s T cas = 1. From the above and equation (8), 
we find 



, T _ (LPp\~* ( X Coul sm 2 6 \-3 _ J_ / Ma \ 4/3 / L \ 2/3 / 10 10 cm \- 2/3 /10 7 K\ 



(A4) 



B. General Nonthermal Damping Rates 

The nonthermal electron damping expression (29) can be integrated for a power law distribution 
of electrons. The general expression is 
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(Bl) 
for a = 3.5, 
for a = 4, 



where 7$ = 1 + Ei/{m e c 2 ) are the Lorentz factors corresponding to kinetic energies Eq and E m . In 
solar case, 7 m -1<1, the above expression then simplifies to equation (31) in the main section. 
Similarly we can get the extreme relativistic limit (7™ S> 1) of equation (32). 
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